Active electromagnetic metamaterial based on spin torque 

oscillators 



Gen Tatara^,^, Hiroaki T. Ueda^, Katsuhisa Taguchi^, 
Yuta Sasaki^, Miyuki Nishijima^, Akihito Takeuchi^ 
^ Graduate School of Science and Engineering, 
Tokyo Metropolitan University, Hachioji, Tokyo 192-0397 Japan 
^ RIKEN Advanced Science Institute, 
2-1 Hirosawa, Wako, Saitama, 351-0198 Japan 
^ Department of Applied Physics, The University of Tokyo, 
7-3-1 Hongo, Tokyo 113-8656 Japan 
(Dated: December 7, 2012) 

Abstract 

We propose theoretically an active material for electromagnetic radiation with frequency of GHz 
by use of spin-torque oscillators. The origin of the amplification is the energy supplied to the 
magnetization by the injected current. We show that close to a resonance with current-driven 
magnetization, the imaginary part of magnetic permeability becomes indeed negative for either of 
the two circular polarizations, resulting in negative imaginary part of refractive index. Besides, the 
real part of the refractive index is also manipulated by the current. Our system thus realizes an 
active filter to obtain circular polarized radiation and/or an electromagnetic metamaterial having 
negative refractive index, both controlled electrically. 



I. INTRODUCTION 



Electric control of material properties is highly important in technologies. For instance, 
manipulation of magnetic structure by electric current is promising for ultra high density 
non- volatile memories [1] and logics [2] , and control of response to electromagnetic radiation 
is useful for sensing, imaging and in communication devices [3]. Here we propose a mech- 
anism of electrically-driven transparency and amplification in ferromagnetic metals in the 
GHz range. The system is a bilayer of thin ferromagnets under a magnetic field and DC 
current. When current is applied to the bilayer, precession of magnetization of the free layer 
is induced [4, 5], resulting in a spin-torque oscillator [6, 7]. The spin-transfer torque induced 
by the applied current acts as a negative damping [5] and then spin-torque oscillator becomes 
an electromagnetic active media close to the resonance frequency. This means that a bilayer 
of ferromagnetic metals, which is a perfect reflector of microwaves, becomes transparent 
when a current is applied. This transparency is a result of active nature of the system, and 
thus the transmitted wave is amplified. The system exhibits in addition a significant effect 
of negative real part of the refractive index, i.e., electromagnetic left-handed metamaterial 
[8]. The present mechanism therefore can be applied to current-induced switch, amplifier, 
polarizer and beamsplitter for GHz waves. 

Negative index of refraction is one of the most counter-intuitive and fascinating phe- 
nomena in the physics of electromagnetic field, where a beam of radiation incident on an 
interface between two materials is refracted towards a wrong direction [9]. Such a possibility 
was theoretically proposed by Veselago many years ago [10]. He discussed that the real part 
of index of refraction becomes negative if the real parts of both permittivity e and magnetic 
permeability /i are negative, and also that negative e and /i arise when a resonance occurs in 
both electric and magnetic properties. In naturally occurring materials, however, negative 
index of refraction is not common. This is because the electric properties of most metallic 
materials are governed by frequencies higher than THz (meV-eV order as the energy scale), 
while frequencies for magnetic properties are lower than lOOGHz, and thus electric and mag- 
netic resonances do not usually occur simultaneously. One way to overcome this difficulty 
is to construct an artificial material, a metamaterial, to realize a system where the plasma 
frequency, of the order of 2400THz (corresponding to 10 eV) in metals, is shifted to lower 
frequencies and/or the magnetic resonance occurs at higher frequencies. A possibility to 
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lower electric resonance was given by Pendry by use of split ring resonators, where induced 
currents in rings result in magnetic resonance [8], and experimental realization of negative 
index of refraction was done in the microwave regime (frequency of 4.8GHz) by combining 
split ring resonators of millimeter size with continuous wires [11]. Later, by use of smaller 
resonators of 20/im size, the frequency was increased to THz regime [12]. Besides split- 
ring resonators, multilayer structure with arrays of holes was used to realize negative index 
of refraction in near- infrared regime (around 150THz) [13]. Very recently, it was demon- 
strated that negative refractive index in the GHz range is realized in a thin film of natural 
ferromagnetic metal [14, 15]. 

For applications, realization of low loss material is of crucial importance. Combining 
negative refraction medium with gain medium (active medium) formed by two-level emitters 
[16] is one possibility [17]. In this paper, we propose an active medium in the GHz regime 
based on a current-driven magnetization of small metallic ferromagnets. The system can 
indeed acquire a negative real part of the refractive index at the same time, realizing an 
active electromagnetic metamaterial with negative refractive index. (In this paper, we call 
systems with negative refractive index the electromagnetic metamaterials.) 

In the case of ferromagnetic resonance driven by an external magnetic field, damping is 
represented by the Gilbert damping constant [18], a, which is positive definite. In other 
words, a static external magnetic field cannot transfer energy to the precessing magnetiza- 
tion, since the motion induced by its torque is always perpendicular to the applied field. 
This situation changes significantly when magnetization dynamics is driven by applying an 
electric current. As pointed out by Slonczewski and Berger [4, 5], when an electric current is 
injected in a layered structure of two ferromagnets, each ferromagnet exerts a torque on the 
magnetization of the other ferromagnet. There occur two types of torques; a spin-transfer 
torque [4, 5] and a field-hke torque [19-22]. The field-like torque acts the same way as a 
magnetic field, while spin-transfer torque is unique for current-induced case. Choosing the 
fixed layer magnetization to be along the precession axis, this torque becomes parallel to 
the damping torque, and thus can reduce the effective damping torque even to a negative 
regime [5]. When the effective damping becomes negative, magnetization starts spontaneous 
precession, resulting in spin-torque oscillators [6] . If an electromagnetic wave with frequency 
close to the precession frequency is injected to such system, the imaginary part of magnetic 
permeability becomes negative. The current-driven magnetization thus works as a gain 
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medium, if the imaginary part is sufficiently large in the negative direction to compensate 
the dissipation in the electric response. We will discuss the real and imaginary parts of the 
refractive index by taking account of the permittivity and present a phase diagram in the 
plane of frequency and applied current. We will demonstrate that a gain medium is indeed 
realized for circularly polarized wave if we tune the frequency close to the resonance by use 
of an external magnetic field. The condition for active material turns out, however, to be 
rather hard to realize; The current density needs to be tuned within the accuracy of 0.1% in 
the worst case. This condition is relaxed by designing a metamaterial to lower the plasma 
frequency like done in Ref. [8]. 

In the next section, we describe a spin-torque oscillator and study its response to an 
incident electromagnetic wave. In Sec. Ill, we briefiy discuss the permittivity, and discuss 
refractive index in Sec. IV. 



II. MAGNETIC PERMEABILITY 



Our system consists of two metallic ferromagnetic layers separated by a thin insulator as 
is a standard setting for current-driven magnetization reversal (Fig. 1(a). One ferromagnet 
is a free layer, whose magnetization (local spins) precesses when current or external field 
is applied, while the other, fixed layer, has a fixed magnetization. When electric current 
is injected perpendicular to the junction, the free layer feels two torques, a spin-transfer 
torque, tst [4, 5], and a field-like torque (perpendicular torque), rg [7, 19, 20]. Denoting 
magnetizations of free and fixed layers by M and Aif, respectively, these torques are 

• 2 • 2 

TST = ^PM X (M X Mf), T^ = ^/3MxMf, (1) 

where j is the applied current density, P is a numerical constant proportional to the spin 
splitting of the conduction electron, ^ is a constant representing the strength of spin relax- 
ation, M = \M\ and a is lattice constant. 

Our first aim is to study the behavior of magnetic permeability when an electromagnetic 
wave is injected into the free layer. We thus include the magnetic field of the electromagnetic 
wave, i?em- The free layer is thinner than the penetration length, which is /im scale for a 
GHz case. For generality, we also apply an external magnetic field, Jfext- In -H'ext, we include 
the effect of magnetic anisotropy field of the system. The equation of motion for M then 
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FIG. 1. (a) Schematic figure of the device consisting of two ferromagnetic layers. Lower layer 
is a ferromagnet with a magnetization M{ fixed in the z direction, and upper layer is a thin 
ferromagnet whose magnetization M can precess when microwave and/or current are applied. 
When the condition for the active media is satisfied for right-handed circular polarization, our 
system works as a active filter to amplify that particular polarization (outgoing green arrow). If 
the system is in the active meta region, i.e., if both real and imaginary parts of refractive index 
are negative, the transmitted wave is circularly polarized and refracted towards a wrong direction 
(red arrow). The system in this regime works as an active beam splitter, (b) A schematic figure of 
torques acting on a local spin with large z component and precessing in the xy plane. The external 
field Hext and fixed layer spin M{ is along z axis, and the magnetic field of the incident microwave, 
i?emi is in the xy plane. The torque due to i^ext, tu{= —jJ-o'yM x iJext), and /3 torque, t/j, point 
towards the tangent to the precession. The damping torque, = -^M x M, tends to tilt the 
precession towards the equilibrium direction along z axis. Spin-transfer torque, tst, acts opposite 
to the damping torque if j < for a right-handed polarized wave, resulting in a negative damping. 



reads [21] 



a 



M = -fiolM X (/fext + -ffcm) + J^M xM+J 



P 

M 



M X {M X Mf) + I3{M X Mi) 



(2) 



where j = ^j, /Uo is the permeabihty in the vacuum, 7(= ^ > 0) is gyromagnetic ratio {g 
is g-factor and electron charge is — e < 0) and a is Gilbert damping constant representing 
relaxation of magnetization. 



We choose the external magnetic field to be along z axis and fixed layer spin as Mf = 
M(0, 0, 1). In this case, the field-like torque points the same direction as the torque from the 
external field, and the dominant part of the spin-transfer torque is parallel or antiparallel to 
the damping torque, as seen from Eq. (2) as follows; iVf x (iVf x iVff) oc iVf x (iVf x i?cxt) oc 
iW X Nl + 0(q;, /3, j, //cm)- The direction of the free magnetization, iVf, is along z direction 
when the electromagnetic field is absent. 

The electromagnetic field is injected parallel to the z axis, namely, has only x and 
y components. When H^xa is applied, the spin thus has a large z component and only small 
components and My in the xy plane, i.e., iVf ~ (M^., M^,, M). We treat M,. and My to the 
linear order. For this geometry, use of H instead of B is convenient, since the components 
of H parallel to the interface are equal for both sides of the interface. We treat the magnetic 
field of the injected electromagnetic field by a plane wave, namely, Jfem oc ^A^-'^-'^^) ^ where 
k is wave vector and uj is angular frequency. Below, we neglect the position dependence, 
considering a small device compared with the wavelength (about 3 mm for a 100 GHz field). 

We consider spin dynamics with angular frequency of w, i.e., M^, oc e"""^*. Defin- 
ing a two-component vector describing spin fluctuation, m = (M^jMy), and hem = 
{Heui,xT Heia,y), the linearized equation reads {ay is a Pauli matrix) (The derivative of 
the z component, M^, is second order in M^, My and i/em-) 

{-co - iPj + (/ioT^^ext + + iuja)ay)m = /j^o^Mayheui- (3) 

The magnetic permeability of the ferromagnct, fi, deflned by fihcm = /Wo^cm + M {fiQ is the 
permeability in the vacuum) is thus obtained as (in a 2 x 2 matrix notation in the xy plane) 



1 + n 



M 



h + — iuja — (a; -|- iPj)ay 



(4) 



where h = no^Hext and Qm = ^^^A*o is an angular frequency due to saturation magnetiza- 
tion. 

We consider a circularly polarized incident wave. The magnetic field h^^ then is /lem = 
/io(l,iO' where Hq is the amplitude and the sign of ± represents left- and right-handed 
polarization (or positive and negative helicity), respectively [23]. Using ayhem — i^iem, we 
obtain the real and imaginary parts of the permeability for each circular polarization, ii±, 
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as 

TuJ + h + /3j 



Re/i± = Ho 



1 + Qm ^ 

{Tco + h + /3j)^ + {ua±Pjy, 



As seen from Eq. (4) and Eq. (5), the resonance arises only for a left-handed light when 
Jfext is along positive z direction. Below, we consider only the left-handed case, 



Re//+ = /lo 



'-UJ + h + iSjy + (ua + Pj) 



Numerical result is shown in Fig. 2. Parameters are chosen as5'=l,a = 2.2A, resulting in 
/XqM ^fXo^S^ 2.2T and Qm = 388GHz, and P = 1 and a = /? = 0.01. 

Let us first look into the imaginary part. We see that when current is zero, it is positive 
definite, indicating that a damping of spin (represented by a) results in a loss of electro- 
magnetic wave. This is no longer true when current is switched on; if Pj < —ua, the 
imaginary part becomes negative for a left-handed wave. The incident wave is therefore 
amplified by coupling to the current-driven magnetization, if the gain overcomes the electric 
loss described by the imaginary part of the permittivity. Remarkably, amplification occurs 
for a particular polarization, either left or right, and is controlled by the current direction, 
as seen from Eq. (5). The present spin-torque oscillator system is thus an active filter to 
obtain a particular circular polarization. The energy necessary for amplification is supplied 
from the applied current. 

Let us look in detail why spin-transfer torque in the present system results in a gain. 
When the fixed magnetization and external field is along z direction, the field-induced torque 
th and current-induced field-like torque rp acting on M are along a tangent to precession 
motion (Fig. 1(b)). The Gilbert damping torque is perpendicular to these torques and tends 
to suppress the precession amplitude, pointing M eventually in the equilibrium direction 
along z axis. The spin-transfer torque in the present configuration is parallel or antiparallel 
to the damping torque depending on the sign of the current. When spin-transfer torque 
points opposite to the damping torque, negative damping is realized resulting in a gain. 

Let us turn to the real part of Eq. (6). As a function of u, the real part of //+ for a 
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FIG. 2. Real and imaginary parts of magnetic permeability for left-handed wave, as function 
of frequency, v = u/{2-k), at /Uo-f^ext = 5T (/i = /xoT-f^ext = 880GHz). Q.m = /"o^ 
The imaginary part becomes negative for a negative current satisfying Eq. (8). 



^ is 388GHz. 



left-handed case has a maximum or minimum of (using a <^1) 



Re/i+(u;^) ~ 1 ± 



(7) 



at w ~ + /3j =p Pj = u^. The real part therefore becomes negative around u = u+ = 
h + f3j + Pj a < Pj < ^ and around u = = h + (3j + \Pj\ if < P j < (see 

Fig. 2(a)). 

Let us look in detail the behavior of /i+ as function of applied current and frequency. 
Equation (6) indicates that negative imaginary part arises when 

~ a 



J < 



(8) 



as depicted by a straight line in Fig. 3(b). It is seen that negative imaginary part, i.e., 
possible gain material, is realized for a broad region of j > 0. This is because the original 
damping constant a is usually small (typically a ^ 0.01). One should note, however, that 
the absolute value of imaginary part is small away from the resonance (Figs. 2 and 3). 
As for the real part of we see from Eq. (6) that it vanishes when 

2 /r^ \ 2 









r 1 2 " 





(PJ? 



n 



M 



(9) 



Here we assumed |/3j| ^ VLm and auj <^ Negative Re/i+ is thus realized inside a oval 

shown in Fig. 3(c). From Eqs. (8) and (9), we see that there are four distinct regions. 
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FIG. 3. Contour plots of (a) real and (b) imaginary parts of /i+ and (c) phase 9^ in a v — j 
plane at ^o-f^ext = 5T. Regions I- IV in (c) are defined as follows. I): Re^+ > and Im/i+ > 0, II): 
Re//+ < 0, Im^+ > 0, III): Re^+ < 0, Im/i+ < and IV): Re/i+ > 0, Imfi^ < 0. The phase is 
defined in the regime —^ + 6e<0^<^+de (see Eq. (20)), where 6^ is a small quantity defined 
by Eq. (16). The phase has a discontinuity between regions III and IV. 

namely, I) Re/i+ > and Im/i+ > 0, II) Refi^ < 0, Imfi^ > 0, III) Refi^ < 0, Im^^ < and 
IV) Re/i+ > 0, Im/i+ < (Fig. 3(c)). 

III. PERMITTIBITY 

So far we discussed solely the magnetic property. Let us now discuss the refractive index, 
n = a/ /ie//ioeo, by including the property of permittivity, e. In the microwave regime, e in 
metals has a large imaginary part due to a strong dissipation and negative real part, reflecting 
the fact that the electromagnetic wave is evanescent. This fact is expressed in equations as 
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follows. The incident electric field E induces a charge current density j — a{u)E, where 
(t{ui) is conductivity, resulting in a decay of E. Combining the Ohm's law with two of the 
Maxwell's equations, V x E = and V x S = /ioj + /ioCoT^, we obtain the permittivity 
modified by the Ohm's law as (using cur <^ 1, where r is the elastic lifetime of electron) 

e = Jl + *^), (10) 

neglecting here the cyclotron motion due to the magnetic field. The imaginary part of e is 
positive, resulting in a finite penetration length of I = (Imy^e^o;)"^ and zero transmission 
of propagating wave [23]. The conductivity in the microwave regime is described by the 
Drude model assuming free electrons, 

a{u) = gb- — ^. — , (11) 

1 — lUlT 

where ctb = — eoUpT is the Boltzmann conductivity and Up = yj^^ is the plasma 
frequency [23]. Including the effect of cyclotron motion, the permittivity is modified to be 
polarization-dependent as [10] (see Appendix A) 

e±^eo(l-^ (12) 

where = ^ — h + Qm the cyclotron angular frequency (electron charge is — e < 0). In 
the low frequency case we consider, a;, Jlc <^ ^, we thus obtain 

Ree± c^eo(^l- (cUpx)^^^^) , (13) 



and 



Im6± c^^^ 60^. (14) 

CO U! 



Therefore, the phase 9^^, defined by e± — |e±|e^^±, is 



^---^-" .[l-(c.,V^- - ^''^ 



Due to strong dissipation by eddy current, 9^ in the GHz range is close to |. Choosing 
kp'^ = O.SA, we have ep = 5.5eV and oup = 1.4 x lO^^Hz {hwp = 9.0eV). Assuming dirty 
metals, we choose epT/h — 10, resulting in uipT — 1.6. For //o-f^xt = 5T and v = u)/{2n) = 



10 



200GHz, we obtain e+/eo ^ 1.0+i x 1.8 x 10^. Since permittivity has large positive imaginary 
part, we define the phase for right-handed polarization, 9^^, as 



The magnitude of deviation, \Se\ is then usually very small; typically of the order of 10~^. 

As seen from Fig. 3(c) and Eq. (9), interesting possibility of negative real part of //+ 
arises in the region ou < flu + h — Qc- Therefore, in this region, the real part of e+ is 
positive and hence 5^ < 0. If the carrier is hole with positive charge, sign of flc reverses, and 
Rce+ < and 6^ > arc realized (if UpT ^1). As we will show in the next section, the sign 
of 6e is important; it determines whether the system is normal (real part of refractive index 
is positive) or meta (negative real part) in the active region. 

IV. REFRACTIVE INDEX 

A. Determination of phase of refractive index 

The Maxwell's equations impose a relation between the wave vector and angular frequency 
as = Cyuw^, where /i and e are permeability and permittivity, respectively. The relative 
refractive index is thus given by 



The sign here, and more precisely the phase of n in the complex plane is crucially important 
in discussing the active nature and sign of Ren. We here demonstrate that the phase is 
determined by the boundary condition imposing that the energy flow is continuous at the 
interface. We consider an interface perpendicular to z axis. Two of the Maxwell's equations, 
V X E = -//^ and V x = then leads to the continuity of E^,Ey, and Hy. We 
consider for simplicity an incident radiation perpendicular to the interface. The first of 
the above two equation reads in terms of wave vector, k = {0,0,k), as k x E = ufiH, 
i.e., {Hx,Hy) = ^{—Ey, Ex), where n = K The energy flow represented by the Poynting 
vector along z axis is — {E x H)^ — '^{E^ + Ey). The sign of the real part oln/ ji thus 
determines the direction of the energy flow, and it needs to be positive to describe the case 
where the incident radiation enters the material. In terms of the phase, defined as e = |e|e'^% 



K = ^ + 



(16) 




(17) 
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/X = |/x|e and n — |n|e this boundary condition requires that 

-l<en-e,<l. (18) 

The case of — — | corresponds to a perfect reflection, since ReP^ vanishes. Since 
Qn— + Oe), the above condition reads 

-n<9^-ee< TT. (19) 

In the present metaUic system, 9^ is close to | in the microwave regime, and thus 9^ = 
I + 6e, where 6^ is a small deviation (positive or negative) as explained in Sec. III. The 
phase 9n is then lies in the region (Fig. 3(c)) 

-| + 5e<^M< Y + (20) 

and 9n takes a value of (Fig. 4(a)) 

S,<9n<7r + 6,. (21) 

From Eq. (21), when > 0, we have three possibilities; 1) normal dissipative, corre- 
sponding to Se < 9n < f ( — I + 5e < < | — 5e), 2) meta dissipative, corresponding to 
I < ^„ < TT ( I — < < ^ — 5e), and 3) meta gain, corresponding to tt < ^„ < tt + 5^ ( 
^ — 5^ < 9^ < ^+6^). Only the third case is the active case, but the window for this region 
is very narrow. When 5^ < 0, wc have 1') normal dissipative, corresponding to < < | 
( — I — Se < 9^ < I — de), 2') meta dissipative, corresponding to | < 6*^ < tt — \Se\ 
( I — 5e < ^/i < ^ + Se), and 3') normal gain, corresponding to —\5e\ < 9n < ( 
— l + ^g < 9/^ < — | — 5g). Again the active case arises in a very narrow window with 
width of A9^ ~ 2S,. 

From the above argument, we see that spin-torque oscillator system realizes negative 
refraction in rather a broad area in the u-j plane, while active material in contrast needs 
fine tuning of applied current to satisfy ^ — 5^ < 9^ < ^ + 5^ (if (5g > 0) or — | — < 
9fj_ < — f + l^el (if < 0). As we have seen in the preceding section, 5^ < in the case 
of electron-like carrier, while 5^ > in the hole-like case, and thus active regime is normal 
(positive Ren) for electron while it is meta (Ren < 0) for hole. 
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FIG. 4. (a) Phase, (b) absolute value, (c) real part and (d) imaginary part of refractive index for 
an incident light with left-handed polarization in a spin-torque oscillator device in the u-j plane. 
Region I and II are Re;U+ > and Im^+ > 0, and Re/i+ < and Im^-|- > 0, respectively. Phase 
is defined in the regime 6^: < 9n < tt + 6^ (Eq. (21)), and there is a discontinuity of phase from 
5^ to IT + 5^ (represented by a thick line in Fig. 4(a)). Active region is realized very close to 
the discontinuity for the region satisfying Eq. (26). (The active region is very narrow and not 
recognizable in the figure). 



B. Refractive index of spin-torque oscillator 



a. Dissipative region Figure 4 shows results of refractive index, n, of spin-torque os- 
cillator. The magnitude of n (Fig. 4(b)) is generally large because of large imaginary part 
of permeability in the GHz range. Negative real part, i.e., electromagnetic metamaterial, is 
realized in the regime inside the oval determined by Eq. (9) (Fig. 4(c)). Most area outside 
of the oval is regime 1) or 1') in the preceding subsection, namely, it is normal dissipative 
(ND) regime. Inside the oval, meta dissipative (MD) is realized corresponding to regime 
2) or 2'). As we have discussed in the preceding subsection, there is a discontinuity of the 

13 



phase Ofj, at the boundary of the oval in the regime j < —au/P from ^ to — |, and active 
material (negative imaginary part) is realized in a very narrow window around the discon- 
tinuity represented by a thick line in Fig. 4(a) (the active regime is not recognizable in the 

figure) . 

Appearance of negative refractive index (electromagnetic meta material) in a broad region 
around j = in Fig. 4(c) is consistent with recent experimental discovery of negative 
refractive index in Co [15]. 

h. Active (gain) region Let us study the active regime (case 3) or 3')) in detail. The 
condition for the active regime is 

K\ ^ (22) 
where 5^ = 6^^^ + |. Since Re/x+ ~ in this regime, 5^ is given by 

From Eq. (9), the current density realizing 5^ = (Re/i+ = 0) for j < is given by 
j = jc = —^[auj + -^Z {uj — /i)(Oc — <^)]- Defining small deviation from jc as J = jc + The 
real and imaginary parts read 

2P{Pl + au) 

Ren+ Ho ~ — — —Sj 

Um{oo - h) 

Im/x+//xo ~ , + 7^ ^—r^Sj, 24 

u — n iim oj — n 



neglecting the higher order of Sj. Thus 



2P 

^ ^Sj. (25) 
The condition for the active media, < \5^\, then becomes 

m < ^2^M\Se\. (26) 

Using the parameters used in Sec. Ill, jp^^M = 6.4 x lO^^A/m^. The window for the active 
regime is thus \6j\ ^ 6 x lO^A/m^ if we use \S^\ = 10~^. Since the applied current density 
is of the order of lO^^A/m^ (see Fig. 4), therefore, fine tuning of current density within 
0.1% is necessary to realize active regime. The active regime is normal (positive Ren) for 
electron-like carrier while it is meta (Ren < 0) for the hole-like case. 
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The high accuracy required for active media is due to strong dissipation in the electric 
response (e) as we have seen. In reahty, this condition is relaxed by several reasons. Most 
obvious one is the lowering of the plasma frequency in dilute systems. In fact, spin-torque 
oscillators are usually made of small thin films with a diameter of lOOnm or less [6]. For 
sufficient electromagnetic response, we need many of such small oscillators put on a substrate, 
resulting in the reduction of the effective electron density in the whole system. The plasma 
frequency, proportional to the square root of the density, is therefore reduced. Reducing the 
effective density by a factor of 0.01 results in enhancement of by a factor of 10, and the 
accuracy needed is relaxed to be 1%. The effective mass of the electric oscillation can be 
tuned, too [8], and this helps the enhancement of \S^\. 

For an experimental realization of the present active material, we need to take into account 
the fact that spin-torque oscillators have thin free layers, with thickness less than lOOnm, 
to excite a coherent oscillation of the magnetization. A single free layer is therefore thinner 
than the penetration depth, which is of the order of micron meter, resulting in almost perfect 
transmission even in the dissipative regime. To observe the active effect we have predicted, 
multilayer structures of spin-torque oscillators where the sum of the free layer is close to fim 
scale is necessary. Besides, to place many spin-torque multilayers as an array is suitable for 
a large signal and also from the viewpoint of reducing dissipation, as we mentioned above. 

V. SUMMARY 

We have shown that spin-torque oscillators act as an active media and/or an electro- 
magnetic metamaterial with negative refractive index for circularly polarized microwaves. 
Without spin-torque oscillators, metals thicker than the penetration depth are almost per- 
fect reflectors of microwaves. What we have shown above is that when the ampliflcation 
by the current-driven magnetization precession in multilayers exceeds the dissipation due to 
the Ohm's law, current-driven ferromagnetic metal becomes transparent in the microwave 
regime. Experimental verification of this current-driven microwave switching is of great 
interest. 

For experimental studies, films with a strong perpendicular easy axis anisotropy like FePt 
and CoPt are of particular interest, since the effective magnetic field (//o-f^ext) due to the 
anisotropy exceeds 10 T in these materials, and thus the present effects are realized at higher 
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frequencies. 



ACKNOWLEDGMENTS 



G.T. thanks S. Tomita, T. Ueda, K. Sawada and M. Gonokami for valuable comments 

and discussions. This work was supported by a Grant-in-Aid for Scientific Research (B) 
(Grant No. 22340104) from Japan Society for the Promotion of Science and UK-Japanese 
GoUaboration on Gurrent- Driven Domain Wall Dynamics from JST. A.T. and K.T. are 
financially supported by the Japan Society for the Promotion of Science for Young Scientists. 



Appendix A: Derivation of permittibity 

Equation of motion of the electron (with charge q = —e) under electric and magnetic 
fields is 

z, = ^ + ^(^xB)-^. (Al) 

mm T 

We consider the case B is along z direction {B — (0, 0, B) and E is in the xy plane. The 
equation then reads 

(1 - iujT)v, ^ ^VyB = ^E,, (A2) 
rn m 

(1 - iujT)vy + —v^B = —Ey, (A3) 
m m 



or 



l-iuT- i^^ay] V = —E (A4) 
m " J m 



For left- and right-haned polarization, E± = -E'o(l, ±i), the velocity is 



qr 



^± = , ■ " .qrB ^, (A5) 

• m 



or 

er 



1 — tur ± 



Using j — —env and P = -j, e = eo + C is therefore obtained as 



1 
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Abstract 

We propose theoretically an active material for electromagnetic radiation with frequency of GHz 
by use of spin-torque oscillators. The origin of the amplification is the energy supplied to the 
magnetization by the injected current. We show that close to a resonance with current-driven 
magnetization, the imaginary part of magnetic permeability becomes indeed negative for either of 
the two circular polarizations, resulting in negative imaginary part of refractive index. Besides, the 
real part of the refractive index is also manipulated by the current. Our system thus realizes an 
active filter to obtain circular polarized radiation and/or an electromagnetic metamaterial having 
negative refractive index, both controlled electrically. 



I. INTRODUCTION 



Electric control of material properties is highly important in technologies. For instance, 
manipulation of magnetic structure by electric current is promising for ultra high density 
non- volatile memories [1] and logics [2] , and control of response to electromagnetic radiation 
is useful for sensing, imaging and in communication devices [3]. Here we propose a mech- 
anism of electrically-driven transparency and amplification in ferromagnetic metals in the 
GHz range. The system is a bilayer of thin ferromagnets under a magnetic field and DC 
current. When current is applied to the bilayer, precession of magnetization of the free layer 
is induced [4, 5], resulting in a spin-torque oscillator [6, 7]. The spin-transfer torque induced 
by the applied current acts as a negative damping [5] and then spin-torque oscillator becomes 
an electromagnetic active media close to the resonance frequency. This means that a bilayer 
of ferromagnetic metals, which is a perfect reflector of microwaves, becomes transparent 
when a current is applied. This transparency is a result of active nature of the system, and 
thus the transmitted wave is amplified. The system exhibits in addition a significant effect 
of negative real part of the refractive index, i.e., electromagnetic left-handed metamaterial 
[8]. The present mechanism therefore can be applied to current-induced switch, amplifier, 
polarizer and beamsplitter for GHz waves. 

Negative index of refraction is one of the most counter-intuitive and fascinating phe- 
nomena in the physics of electromagnetic field, where a beam of radiation incident on an 
interface between two materials is refracted towards a wrong direction [9]. Such a possibility 
was theoretically proposed by Veselago many years ago [10]. He discussed that the real part 
of index of refraction becomes negative if the real parts of both permittivity e and magnetic 
permeability /i are negative, and also that negative e and /i arise when a resonance occurs in 
both electric and magnetic properties. In naturally occurring materials, however, negative 
index of refraction is not common. This is because the electric properties of most metallic 
materials are governed by frequencies higher than THz (meV-eV order as the energy scale), 
while frequencies for magnetic properties are lower than lOOGHz, and thus electric and mag- 
netic resonances do not usually occur simultaneously. One way to overcome this difficulty 
is to construct an artificial material, a metamaterial, to realize a system where the plasma 
frequency, of the order of 2400THz (corresponding to 10 eV) in metals, is shifted to lower 
frequencies and/or the magnetic resonance occurs at higher frequencies. A possibility to 
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lower electric resonance was given by Pendry by use of split ring resonators, where induced 
currents in rings result in magnetic resonance [8], and experimental realization of negative 
index of refraction was done in the microwave regime (frequency of 4.8GHz) by combining 
split ring resonators of millimeter size with continuous wires [11]. Later, by use of smaller 
resonators of 20/im size, the frequency was increased to THz regime [12]. Besides split- 
ring resonators, multilayer structure with arrays of holes was used to realize negative index 
of refraction in near- infrared regime (around 150THz) [13]. Very recently, it was demon- 
strated that negative refractive index in the GHz range is realized in a thin film of natural 
ferromagnetic metal [14, 15]. 

For applications, realization of low loss material is of crucial importance. Combining 
negative refraction medium with gain medium (active medium) formed by two-level emitters 
[16] is one possibility [17]. In this paper, we propose an active medium in the GHz regime 
based on a current-driven magnetization of small metallic ferromagnets. The system can 
indeed acquire a negative real part of the refractive index at the same time, realizing an 
active electromagnetic metamaterial with negative refractive index. (In this paper, we call 
systems with negative refractive index the electromagnetic metamaterials.) 

In the case of ferromagnetic resonance driven by an external magnetic field, damping is 
represented by the Gilbert damping constant [18], a, which is positive definite. In other 
words, a static external magnetic field cannot transfer energy to the precessing magnetiza- 
tion, since the motion induced by its torque is always perpendicular to the applied field. 
This situation changes significantly when magnetization dynamics is driven by applying an 
electric current. As pointed out by Slonczewski and Berger [4, 5], when an electric current is 
injected in a layered structure of two ferromagnets, each ferromagnet exerts a torque on the 
magnetization of the other ferromagnet. There occur two types of torques; a spin-transfer 
torque [4, 5] and a field-hke torque [19-22]. The field-like torque acts the same way as a 
magnetic field, while spin-transfer torque is unique for current-induced case. Choosing the 
fixed layer magnetization to be along the precession axis, this torque becomes parallel to 
the damping torque, and thus can reduce the effective damping torque even to a negative 
regime [5]. When the effective damping becomes negative, magnetization starts spontaneous 
precession, resulting in spin-torque oscillators [6] . If an electromagnetic wave with frequency 
close to the precession frequency is injected to such system, the imaginary part of magnetic 
permeability becomes negative. The current-driven magnetization thus works as a gain 
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medium, if the imaginary part is sufficiently large in the negative direction to compensate 
the dissipation in the electric response. We will discuss the real and imaginary parts of the 
refractive index by taking account of the permittivity and present a phase diagram in the 
plane of frequency and applied current. We will demonstrate that a gain medium is indeed 
realized for circularly polarized wave if we tune the frequency close to the resonance by use 
of an external magnetic field. The condition for active material turns out, however, to be 
rather hard to realize; The current density needs to be tuned within the accuracy of 0.1% in 
the worst case. This condition is relaxed by designing a metamaterial to lower the plasma 
frequency like done in Ref. [8]. 

In the next section, we describe a spin-torque oscillator and study its response to an 
incident electromagnetic wave. In Sec. Ill, we briefiy discuss the permittivity, and discuss 
refractive index in Sec. IV. 



II. MAGNETIC PERMEABILITY 



Our system consists of two metallic ferromagnetic layers separated by a thin insulator as 
is a standard setting for current-driven magnetization reversal (Fig. 1(a). One ferromagnet 
is a free layer, whose magnetization (local spins) precesses when current or external field 
is applied, while the other, fixed layer, has a fixed magnetization. When electric current 
is injected perpendicular to the junction, the free layer feels two torques, a spin-transfer 
torque, tst [4, 5], and a field-like torque (perpendicular torque), rg [7, 19, 20]. Denoting 
magnetizations of free and fixed layers by M and Aif, respectively, these torques are 

• 2 • 2 

TST = ^PM X (M X Mf), T^ = ^/3MxMf, (1) 

where j is the applied current density, P is a numerical constant proportional to the spin 
splitting of the conduction electron, ^ is a constant representing the strength of spin relax- 
ation, M = \M\ and a is lattice constant. 

Our first aim is to study the behavior of magnetic permeability when an electromagnetic 
wave is injected into the free layer. We thus include the magnetic field of the electromagnetic 
wave, i?em- The free layer is thinner than the penetration length, which is /im scale for a 
GHz case. For generality, we also apply an external magnetic field, Jfext- In -H'ext, we include 
the effect of magnetic anisotropy field of the system. The equation of motion for M then 
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FIG. 1. (a) Schematic figure of the device consisting of two ferromagnetic layers. Lower layer 
is a ferromagnet with a magnetization M{ fixed in the z direction, and upper layer is a thin 
ferromagnet whose magnetization M can precess when microwave and/or current are applied. 
When the condition for the active media is satisfied for right-handed circular polarization, our 
system works as a active filter to amplify that particular polarization (outgoing green arrow). If 
the system is in the active meta region, i.e., if both real and imaginary parts of refractive index 
are negative, the transmitted wave is circularly polarized and refracted towards a wrong direction 
(red arrow). The system in this regime works as an active beam splitter, (b) A schematic figure of 
torques acting on a local spin with large z component and precessing in the xy plane. The external 
field Hext and fixed layer spin M{ is along z axis, and the magnetic field of the incident microwave, 
i?emi is in the xy plane. The torque due to i^ext, tu{= —jJ-o'yM x iJext), and /3 torque, t/j, point 
towards the tangent to the precession. The damping torque, = -^M x M, tends to tilt the 
precession towards the equilibrium direction along z axis. Spin-transfer torque, tst, acts opposite 
to the damping torque if j < for a right-handed polarized wave, resulting in a negative damping. 



reads [21] 



a 



M = -fiolM X (/fext + -ffcm) + J^M xM+J 



P 

M 



M X {M X Mf) + I3{M X Mi) 



(2) 



where j = ^j, /Uo is the permeabihty in the vacuum, 7(= ^ > 0) is gyromagnetic ratio {g 
is g-factor and electron charge is — e < 0) and a is Gilbert damping constant representing 
relaxation of magnetization. 



We choose the external magnetic field to be along z axis and fixed layer spin as Mf = 
M(0, 0, 1). In this case, the field-like torque points the same direction as the torque from the 
external field, and the dominant part of the spin-transfer torque is parallel or antiparallel to 
the damping torque, as seen from Eq. (2) as follows; iVf x (iVf x iVff) oc iVf x (iVf x i?cxt) oc 
iW X Nl + 0(q;, /3, j, //cm)- The direction of the free magnetization, iVf, is along z direction 
when the electromagnetic field is absent. 

The electromagnetic field is injected parallel to the z axis, namely, has only x and 
y components. When H^xa is applied, the spin thus has a large z component and only small 
components and My in the xy plane, i.e., iVf ~ (M^., M^,, M). We treat M,. and My to the 
linear order. For this geometry, use of H instead of B is convenient, since the components 
of H parallel to the interface are equal for both sides of the interface. We treat the magnetic 
field of the injected electromagnetic field by a plane wave, namely, Jfem oc ^A^-'^-'^^) ^ where 
k is wave vector and uj is angular frequency. Below, we neglect the position dependence, 
considering a small device compared with the wavelength (about 3 mm for a 100 GHz field). 

We consider spin dynamics with angular frequency of w, i.e., M^, oc e"""^*. Defin- 
ing a two-component vector describing spin fluctuation, m = (M^jMy), and hem = 
{Heui,xT Heia,y), the linearized equation reads {ay is a Pauli matrix) (The derivative of 
the z component, M^, is second order in M^, My and i/em-) 

{-co - iPj + (/ioT^^ext + + iuja)ay)m = /j^o^Mayheui- (3) 

The magnetic permeability of the ferromagnct, fi, deflned by fihcm = /Wo^cm + M {fiQ is the 
permeability in the vacuum) is thus obtained as (in a 2 x 2 matrix notation in the xy plane) 



1 + n 



M 



h + — iuja — (a; -|- iPj)ay 



(4) 



where h = no^Hext and Qm = ^^^A*o is an angular frequency due to saturation magnetiza- 
tion. 

We consider a circularly polarized incident wave. The magnetic field h^^ then is /lem = 
/io(l,iO' where Hq is the amplitude and the sign of ± represents left- and right-handed 
polarization (or positive and negative helicity), respectively [23]. Using ayhem — i^iem, we 
obtain the real and imaginary parts of the permeability for each circular polarization, ii±, 
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as 

TuJ + h + /3j 



Re/i± = Ho 



1 + Qm ^ 

{Tco + h + /3j)^ + {ua±Pjy, 



As seen from Eq. (4) and Eq. (5), the resonance arises only for a left-handed light when 
Jfext is along positive z direction. Below, we consider only the left-handed case, 



Re//+ = /lo 



'-UJ + h + iSjy + (ua + Pj) 



Numerical result is shown in Fig. 2. Parameters are chosen as5'=l,a = 2.2A, resulting in 
/XqM ^fXo^S^ 2.2T and Qm = 388GHz, and P = 1 and a = /? = 0.01. 

Let us first look into the imaginary part. We see that when current is zero, it is positive 
definite, indicating that a damping of spin (represented by a) results in a loss of electro- 
magnetic wave. This is no longer true when current is switched on; if Pj < —ua, the 
imaginary part becomes negative for a left-handed wave. The incident wave is therefore 
amplified by coupling to the current-driven magnetization, if the gain overcomes the electric 
loss described by the imaginary part of the permittivity. Remarkably, amplification occurs 
for a particular polarization, either left or right, and is controlled by the current direction, 
as seen from Eq. (5). The present spin-torque oscillator system is thus an active filter to 
obtain a particular circular polarization. The energy necessary for amplification is supplied 
from the applied current. 

Let us look in detail why spin-transfer torque in the present system results in a gain. 
When the fixed magnetization and external field is along z direction, the field-induced torque 
th and current-induced field-like torque rp acting on M are along a tangent to precession 
motion (Fig. 1(b)). The Gilbert damping torque is perpendicular to these torques and tends 
to suppress the precession amplitude, pointing M eventually in the equilibrium direction 
along z axis. The spin-transfer torque in the present configuration is parallel or antiparallel 
to the damping torque depending on the sign of the current. When spin-transfer torque 
points opposite to the damping torque, negative damping is realized resulting in a gain. 

Let us turn to the real part of Eq. (6). As a function of u, the real part of //+ for a 
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FIG. 2. Real and imaginary parts of magnetic permeability for left-handed wave, as function 
of frequency, v = u/{2-k), at /Uo-f^ext = 5T (/i = /xoT-f^ext = 880GHz). Q.m = /"o^ 
The imaginary part becomes negative for a negative current satisfying Eq. (8). 



^ is 388GHz. 



left-handed case has a maximum or minimum of (using a <^1) 



Re/i+(u;^) ~ 1 ± 



(7) 



at w ~ + /3j =p Pj = u^. The real part therefore becomes negative around u = u+ = 
h + f3j + Pj a < Pj < ^ and around u = = h + (3j + \Pj\ if < P j < (see 

Fig. 2(a)). 

Let us look in detail the behavior of /i+ as function of applied current and frequency. 
Equation (6) indicates that negative imaginary part arises when 

~ a 



J < 



(8) 



as depicted by a straight line in Fig. 3(b). It is seen that negative imaginary part, i.e., 
possible gain material, is realized for a broad region of j > 0. This is because the original 
damping constant a is usually small (typically a ^ 0.01). One should note, however, that 
the absolute value of imaginary part is small away from the resonance (Figs. 2 and 3). 
As for the real part of we see from Eq. (6) that it vanishes when 

2 /r^ \ 2 









r 1 2 " 





(PJ? 



n 



M 



(9) 



Here we assumed |/3j| ^ VLm and auj <^ Negative Re/i+ is thus realized inside a oval 

shown in Fig. 3(c). From Eqs. (8) and (9), we see that there are four distinct regions. 
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FIG. 3. Contour plots of (a) real and (b) imaginary parts of /i+ and (c) phase 9^ in a v — j 
plane at ^o-f^ext = 5T. Regions I- IV in (c) are defined as follows. I): Re^+ > and Im/i+ > 0, II): 
Re//+ < 0, Im^+ > 0, III): Re^+ < 0, Im/i+ < and IV): Re/i+ > 0, Imfi^ < 0. The phase is 
defined in the regime —^ + 6e<0^<^+de (see Eq. (20)), where 6^ is a small quantity defined 
by Eq. (16). The phase has a discontinuity between regions III and IV. 

namely, I) Re/i+ > and Im/i+ > 0, II) Refi^ < 0, Imfi^ > 0, III) Refi^ < 0, Im^^ < and 
IV) Re/i+ > 0, Im/i+ < (Fig. 3(c)). 

III. PERMITTIBITY 

So far we discussed solely the magnetic property. Let us now discuss the refractive index, 
n = a/ /ie//ioeo, by including the property of permittivity, e. In the microwave regime, e in 
metals has a large imaginary part due to a strong dissipation and negative real part, reflecting 
the fact that the electromagnetic wave is evanescent. This fact is expressed in equations as 
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follows. The incident electric field E induces a charge current density j — a{u)E, where 
(t{ui) is conductivity, resulting in a decay of E. Combining the Ohm's law with two of the 
Maxwell's equations, V x E = and V x S = /ioj + /ioCoT^, we obtain the permittivity 
modified by the Ohm's law as (using cur <^ 1, where r is the elastic lifetime of electron) 

e = Jl + *^), (10) 

neglecting here the cyclotron motion due to the magnetic field. The imaginary part of e is 
positive, resulting in a finite penetration length of I = (Imy^e^o;)"^ and zero transmission 
of propagating wave [23]. The conductivity in the microwave regime is described by the 
Drude model assuming free electrons, 

a{u) = gb- — ^. — , (11) 

1 — lUlT 

where ctb = — eoUpT is the Boltzmann conductivity and Up = yj^^ is the plasma 
frequency [23]. Including the effect of cyclotron motion, the permittivity is modified to be 
polarization-dependent as [10] (see Appendix A) 

e±^eo(l-^ (12) 

where = ^ — h + Qm the cyclotron angular frequency (electron charge is — e < 0). In 
the low frequency case we consider, a;, Jlc <^ ^, we thus obtain 

Ree± c^eo(^l- (cUpx)^^^^) , (13) 



and 



Im6± c^^^ 60^. (14) 

CO U! 



Therefore, the phase 9^^, defined by e± — |e±|e^^±, is 



^---^-" .[l-(c.,V^- - ^''^ 



Due to strong dissipation by eddy current, 9^ in the GHz range is close to |. Choosing 
kp'^ = O.SA, we have ep = 5.5eV and oup = 1.4 x lO^^Hz {hwp = 9.0eV). Assuming dirty 
metals, we choose epT/h — 10, resulting in uipT — 1.6. For //o-f^xt = 5T and v = u)/{2n) = 
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200GHz, we obtain e+/eo ^ 1.0+i x 1.8 x 10^. Since permittivity has large positive imaginary 
part, we define the phase for right-handed polarization, 9^^, as 



The magnitude of deviation, \Se\ is then usually very small; typically of the order of 10~^. 

As seen from Fig. 3(c) and Eq. (9), interesting possibility of negative real part of //+ 
arises in the region ou < flu + h — Qc- Therefore, in this region, the real part of e+ is 
positive and hence 5^ < 0. If the carrier is hole with positive charge, sign of flc reverses, and 
Rce+ < and 6^ > arc realized (if UpT ^1). As we will show in the next section, the sign 
of 6e is important; it determines whether the system is normal (real part of refractive index 
is positive) or meta (negative real part) in the active region. 

IV. REFRACTIVE INDEX 

A. Determination of phase of refractive index 

The Maxwell's equations impose a relation between the wave vector and angular frequency 
as = Cyuw^, where /i and e are permeability and permittivity, respectively. The relative 
refractive index is thus given by 



The sign here, and more precisely the phase of n in the complex plane is crucially important 
in discussing the active nature and sign of Ren. We here demonstrate that the phase is 
determined by the boundary condition imposing that the energy flow is continuous at the 
interface. We consider an interface perpendicular to z axis. Two of the Maxwell's equations, 
V X E = -//^ and V x = then leads to the continuity of E^,Ey, and Hy. We 
consider for simplicity an incident radiation perpendicular to the interface. The first of 
the above two equation reads in terms of wave vector, k = {0,0,k), as k x E = ufiH, 
i.e., {Hx,Hy) = ^{—Ey, Ex), where n = K The energy flow represented by the Poynting 
vector along z axis is — {E x H)^ — '^{E^ + Ey). The sign of the real part oln/ ji thus 
determines the direction of the energy flow, and it needs to be positive to describe the case 
where the incident radiation enters the material. In terms of the phase, defined as e = |e|e'^% 



K = ^ + 



(16) 




(17) 
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/X = |/x|e and n — |n|e this boundary condition requires that 

-l<en-e,<l. (18) 

The case of — — | corresponds to a perfect reflection, since ReP^ vanishes. Since 
Qn— + Oe), the above condition reads 

-n<9^-ee< TT. (19) 

In the present metaUic system, 9^ is close to | in the microwave regime, and thus 9^ = 
I + 6e, where 6^ is a small deviation (positive or negative) as explained in Sec. III. The 
phase 9n is then lies in the region (Fig. 3(c)) 

-| + 5e<^M< Y + (20) 

and 9n takes a value of (Fig. 4(a)) 

S,<9n<7r + 6,. (21) 

From Eq. (21), when > 0, we have three possibilities; 1) normal dissipative, corre- 
sponding to Se < 9n < f ( — I + 5e < < | — 5e), 2) meta dissipative, corresponding to 
I < ^„ < TT ( I — < < ^ — 5e), and 3) meta gain, corresponding to tt < ^„ < tt + 5^ ( 
^ — 5^ < 9^ < ^+6^). Only the third case is the active case, but the window for this region 
is very narrow. When 5^ < 0, wc have 1') normal dissipative, corresponding to < < | 
( — I — Se < 9^ < I — de), 2') meta dissipative, corresponding to | < 6*^ < tt — \Se\ 
( I — 5e < ^/i < ^ + Se), and 3') normal gain, corresponding to —\5e\ < 9n < ( 
— l + ^g < 9/^ < — | — 5g). Again the active case arises in a very narrow window with 
width of A9^ ~ 2S,. 

From the above argument, we see that spin-torque oscillator system realizes negative 
refraction in rather a broad area in the u-j plane, while active material in contrast needs 
fine tuning of applied current to satisfy ^ — 5^ < 9^ < ^ + 5^ (if (5g > 0) or — | — < 
9fj_ < — f + l^el (if < 0). As we have seen in the preceding section, 5^ < in the case 
of electron-like carrier, while 5^ > in the hole-like case, and thus active regime is normal 
(positive Ren) for electron while it is meta (Ren < 0) for hole. 
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FIG. 4. (a) Phase, (b) absolute value, (c) real part and (d) imaginary part of refractive index for 
an incident light with left-handed polarization in a spin-torque oscillator device in the u-j plane. 
Region I and II are Re;U+ > and Im^+ > 0, and Re/i+ < and Im^-|- > 0, respectively. Phase 
is defined in the regime 6^: < 9n < tt + 6^ (Eq. (21)), and there is a discontinuity of phase from 
5^ to IT + 5^ (represented by a thick line in Fig. 4(a)). Active region is realized very close to 
the discontinuity for the region satisfying Eq. (26). (The active region is very narrow and not 
recognizable in the figure). 



B. Refractive index of spin-torque oscillator 



a. Dissipative region Figure 4 shows results of refractive index, n, of spin-torque os- 
cillator. The magnitude of n (Fig. 4(b)) is generally large because of large imaginary part 
of permeability in the GHz range. Negative real part, i.e., electromagnetic metamaterial, is 
realized in the regime inside the oval determined by Eq. (9) (Fig. 4(c)). Most area outside 
of the oval is regime 1) or 1') in the preceding subsection, namely, it is normal dissipative 
(ND) regime. Inside the oval, meta dissipative (MD) is realized corresponding to regime 
2) or 2'). As we have discussed in the preceding subsection, there is a discontinuity of the 
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phase Ofj, at the boundary of the oval in the regime j < —au/P from ^ to — |, and active 
material (negative imaginary part) is realized in a very narrow window around the discon- 
tinuity represented by a thick line in Fig. 4(a) (the active regime is not recognizable in the 

figure) . 

Appearance of negative refractive index (electromagnetic meta material) in a broad region 
around j = in Fig. 4(c) is consistent with recent experimental discovery of negative 
refractive index in Co [15]. 

h. Active (gain) region Let us study the active regime (case 3) or 3')) in detail. The 
condition for the active regime is 

K\ ^ (22) 
where 5^ = 6^^^ + |. Since Re/x+ ~ in this regime, 5^ is given by 

From Eq. (9), the current density realizing 5^ = (Re/i+ = 0) for j < is given by 
j = jc = —^[auj + -^Z {uj — /i)(Oc — <^)]- Defining small deviation from jc as J = jc + The 
real and imaginary parts read 

2P{Pl + au) 

Ren+ Ho ~ — — —Sj 

Um{oo - h) 

Im/x+//xo ~ , + 7^ ^—r^Sj, 24 

u — n iim oj — n 



neglecting the higher order of Sj. Thus 



2P 

^ ^Sj. (25) 
The condition for the active media, < \5^\, then becomes 

m < ^2^M\Se\. (26) 

Using the parameters used in Sec. Ill, jp^^M = 6.4 x lO^^A/m^. The window for the active 
regime is thus \6j\ ^ 6 x lO^A/m^ if we use \S^\ = 10~^. Since the applied current density 
is of the order of lO^^A/m^ (see Fig. 4), therefore, fine tuning of current density within 
0.1% is necessary to realize active regime. The active regime is normal (positive Ren) for 
electron-like carrier while it is meta (Ren < 0) for the hole-like case. 
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The high accuracy required for active media is due to strong dissipation in the electric 
response (e) as we have seen. In reahty, this condition is relaxed by several reasons. Most 
obvious one is the lowering of the plasma frequency in dilute systems. In fact, spin-torque 
oscillators are usually made of small thin films with a diameter of lOOnm or less [6]. For 
sufficient electromagnetic response, we need many of such small oscillators put on a substrate, 
resulting in the reduction of the effective electron density in the whole system. The plasma 
frequency, proportional to the square root of the density, is therefore reduced. Reducing the 
effective density by a factor of 0.01 results in enhancement of by a factor of 10, and the 
accuracy needed is relaxed to be 1%. The effective mass of the electric oscillation can be 
tuned, too [8], and this helps the enhancement of \S^\. 

For an experimental realization of the present active material, we need to take into account 
the fact that spin-torque oscillators have thin free layers, with thickness less than lOOnm, 
to excite a coherent oscillation of the magnetization. A single free layer is therefore thinner 
than the penetration depth, which is of the order of micron meter, resulting in almost perfect 
transmission even in the dissipative regime. To observe the active effect we have predicted, 
multilayer structures of spin-torque oscillators where the sum of the free layer is close to fim 
scale is necessary. Besides, to place many spin-torque multilayers as an array is suitable for 
a large signal and also from the viewpoint of reducing dissipation, as we mentioned above. 

V. SUMMARY 

We have shown that spin-torque oscillators act as an active media and/or an electro- 
magnetic metamaterial with negative refractive index for circularly polarized microwaves. 
Without spin-torque oscillators, metals thicker than the penetration depth are almost per- 
fect reflectors of microwaves. What we have shown above is that when the ampliflcation 
by the current-driven magnetization precession in multilayers exceeds the dissipation due to 
the Ohm's law, current-driven ferromagnetic metal becomes transparent in the microwave 
regime. Experimental verification of this current-driven microwave switching is of great 
interest. 

For experimental studies, films with a strong perpendicular easy axis anisotropy like FePt 
and CoPt are of particular interest, since the effective magnetic field (//o-f^ext) due to the 
anisotropy exceeds 10 T in these materials, and thus the present effects are realized at higher 
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frequencies. 
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Appendix A: Derivation of permittibity 

Equation of motion of the electron (with charge q = —e) under electric and magnetic 
fields is 

z, = ^ + ^(^xB)-^. (Al) 

mm T 

We consider the case B is along z direction {B — (0, 0, B) and E is in the xy plane. The 
equation then reads 

(1 - iujT)v, ^ ^VyB = ^E,, (A2) 
rn m 

(1 - iujT)vy + —v^B = —Ey, (A3) 
m m 



or 



l-iuT- i^^ay] V = —E (A4) 
m " J m 



For left- and right-haned polarization, E± = -E'o(l, ±i), the velocity is 



qr 



^± = , ■ " .qrB ^, (A5) 

• m 



or 

er 



1 — tur ± 



Using j — —env and P = -j, e = eo + C is therefore obtained as 
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